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ABSTRACT 

We propose a new method of analysing a variable component for type 1 active galactic nuclei (AGNs) 
in the near-infrared wavelength region. This analysis uses a multiple regression technique and divides 
the variable component into two components originating in the accretion disk at the center of AGNs 
and from the dust torus that far surrounds the disk. Applying this analysis to the long-term VHK 
monitoring data of MCG+08-11-011 that were obtained by the MAGNUM project, we found that 
the (H — iQ-color temperature of the dust component is T — 1635K ±20K, which agrees with the 
sublimation temperature of dust grains, and that the time delay of K to H variations is At w 6 
days, which indicates the existence of a radial temperature gradient in the dust torus. As for the 
disk component, we found that the power-law spectrum of f u oc v a in the V to near-infrared HK 

bands varies with a fixed index of a « —0.1 hO.4, which is broadly consistent with the irradiated 

standard disk model. The outer part of the disk therefore extends out to a radial distance where the 
temperature decreases to radiate the light in the near-infrared. 

Subject headings: galaxies: active — galaxies: individual (MCG+08-11-011) — accretion disks — dust 



1. INTRODUCTION 

The existence of massive black holes at the center of 
accretion disks in AGNs is well accepted in current AGN 
research. The so-called big blue bump, which is observed 
in the UV to optical wavelength region of spectral energy 
distribution (SED) of many AGNs, is regarded as an evi- 
dence of the accretion disk (e.g., Shields 1978; Malkan & 
Sargent 1982; Czerny & Elvis 1987). Based on another 
distinct feature of near-infrared bumps in their SED, the 
existence of a surrounding dust torus, which is heated 
to the sublimation temperature of dust by the accretion 
disk, is also accepted now (e.g., Sanders et al. 1989; 
Barvainis 1987; Kobayashi et al. 1993). 

In practice, however, it is known that composite SEDs 
constructed from a huge number of AGNs obtained by 
large surveys show their optical colors to be significantly 
redder than predicted by the standard accretion disk 
model (see Francis et al. 1991 for the composite SED 
from the low-redshift bright quasar survey (LBQS); Van- 
den Berk et al. 2001 for the composite SED from the 
Sloan digital sky survey (SDSS)). Since AGNs are super- 
imposed on host galaxies, how serious this discrepancy 
is depends on how reliably the nuclear SED has been 
separated from the SED of the galaxy component. 

In addition to the problem of SED decomposition, a 
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possible contamination of near-infrared light from accre- 
tion disk and dust torus further complicates the inter- 
pretation of the SED. In fact, it is reported that the 
SED of the accretion disk extends to the near-infrared 
wavelength region (Kishimoto et al. 2005; Minezaki et 
al. 2006), and that the temperature gradient exists in 
the dust torus of which the inner part is hotter than the 
outer part (e.g. Barvainis 1987; Glass 1992; Alloin et 
al. 1995). These two factors cause similar effects such 
that variations at shorter near-infrared wavelengths are 
less delayed behind the optical variations, when com- 
pared to variations at longer near-infrared wavelengths. 
Therefore, unless the above factors are taken into ac- 
count simultaneously in the analysis, the time delay of 
near-infrared variations behind the optical is liable to be 
underestimated. 

In order to overcome such complications as well as 
the problem of SED decomposition, we propose a new 
method of analysing the variable SEDs of accretion disks 
and dust tori, based on a multiple regression technique. 
In section 2, we describe the assumptions and methods 
of the analysis. In section 3, we apply the analysis to 
the long-term VHK monitoring data of MCG+08-11- 
011 obtained by the MAGNUM project (Yoshii 2002; 
Kobayashi et al. 2003; Yoshii et al. 2003). Finally, 
in section 4, we describe the applicability of the analysis 
and discuss the results. 

2. MULTIPLE REGRESSION ANALYSIS 

2.1. Overview of the Analysis 

Let F\(t) be the observed flux at wavelength A and 
time t, then it is decomposed as 

F\(t) = F v x ar ' dlsk (t) + F^ ar ' dust (t) + F A const , (1) 

where the first two are the variable components from 
the accretion disk and dust torus, respectively, and the 
last component is a constant, partly offset including the 
contribution from the host galaxy. 
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Kobayashi et al. (1993) showed that all 14 QSOs in 
their sample have almost the same blackbody spectra of 
1500K. This indicates that the flux in the near-infrared 
wavelength region is mainly radiated from the innermost 
part of the dust torus, where the temperature is near the 
sublimation temperature of dust grains. Hence, as we 
observe for each AGN, the near-infrared light curves in 
different bands have almost the same shape and overlap 
each other by shifting the lag time between the variations 
in different bands (e.g., Suganuma et al. 2006). When 
there is a short lag time of r between variations in two 
near-infrared bands, H (1.1 fim) and K(2.2fim) for exam- 
ple, which are intrinsic to the dust torus, the situation 
considered here corresponds to 

F H (t- t)/F k (t) = a. (2) 

Winkler et al. (1992) and Winkler (1997) claimed that 
the optical colors of the variable component for Seyfert 
galaxies remain constant as their brightness changes, 
while recent studies report that QSOs get bluer as they 
get brighter (e.g., Trevese & Vagnetti 2002; Giveon et 
al. 1999). For our analysis focusing on Seyfert galax- 
ies below, we then assume that the variable component 
synchronized with optical variation, which comes from 
the accretion disk, has constant optical V(0.55/J,m) to 
near- infrared H and K colors, i.e., 



F v H ar ' disk (t)/F v {t) 
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\t)/F v (t) = c. 

(3) 

Substitution of equations 2 and 3 in equation 1 therefore 
gives 

F H (t-T)-bxF v (t-T) = ax(F K (t)-cxF v (t))+d, (4) 

where a, b, c, and d are the constant parameters to be de- 
termined by the fit to the VHK monitoring data, based 
on a multiple regression technique described below. 

2.2. Details of Calculation 

For calculation of the multiple regression, Jefferys 
(1980) provides the most generalized multivariate least- 
square process. This process allows errors to appear in 
all variables, which is desirable for astronomical observa- 
tions. In practice, we adopted the Mardquardt method 
described in Jefferys (1981), which is robust for the case 
that the initial trial for the fit is far from the true solu- 
tion. 

The calculation of the regression with a given value 
of lag time r requires interpolation of the monitoring 
data in time, because such data are usually discrete. We 
here adopt the linear interpolation scheme, and smaller 
weights are assigned to the data produced with large 
interpolation. We estimate the uncertainties of the in- 
terpolation, using the structure function (SF) of AGN 
variabilities (Collier & Peterson 2001; Suganuma et al. 
2004): 

SF (^M^ l/( ' ,)_/te)f_y ' (5) 



2 

CT SF 



where the summation runs over all -pairs, I =tj — U 
is the length of interpolation, N(l) is the number of pairs 
at I, f(t) is the flux at time t, and a is the mean observa- 
tional error for the data used for the SF calculation. We 
adopted a length nearer to the neighboring data point 



as the length of interpolation I. While we naturally use 
the observational error as the weight for the data with 
no interpolation, we adopt the larger of the SF error or 
the observational error as representing the uncertainty 
of the interpolated data in the regression analysis. In 
order to determine whether or not the fit well repre- 
sents the observations, the x 2 -tcst is performed. The 
X 2 -value is computed by doubling equation 8 of Jefferys 
(1980), and this doubled value follows the \ 2 distribution 
of N — m — 4 degrees of freedom, where m is the number 
of data sets. 

Finally, the value of lag time r that gives the best fit is 
estimated at the peak of multiple correlation coefficient 
that measures the strength of association between the 
independent and dependent variables. 
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Fig. 1. — Light curves of MCG+08-11-011 for the period from 
8 September 2001 to 31 December 2002 (UTC). Shown arc the 
monitoring data for the B band (open circles), V band (crosses), 
H band (open rectangles), and K band (filled rectangles), where 
the flux from the host galaxy has been subtracted. 



3. APPLICATION OF THE ANALYSIS 

The BVHK monitoring data for MCG+08-11-011 
were obtained by the MAGNUM project in the period 
from 8 September 2001 to 31 December 2002 (Figure 1). 
Aperture photometry was used throughout. An aper- 
ture size of 8.3 arcsec was adopted for all images. The 
flux of MCG+08-11-011 was measured relative to a ref- 
erence star (a 20 oo = 05 h 54 m 47. s 9, (5 2 ooo = +46°23'5") 
whose flux was absolutely calibrated using Landolt stan- 
dard stars (Landolt 1992) in BV and Hunt standard stars 
(Hunt et al. 1997) in HK during the monitoring ob- 
servations. The host galaxy flux was estimated using 
the GALFIT program (Peng et al. 2002) and was sub- 
tracted from the measured flux of MCG+08-11-011. The 
estimated magnitudes of the host galaxy were 16.5 mag 
at B, 15.5 mag at V, 12.0 mag at H and 11.5 mag at K. 
A small error of 0.003 mag caused by the seeing effect, 
arising from the profile difference between the AGN and 
the reference star, was taken into account. Galactic ex- 
tinction was corrected according to Schlegcl, Finkbciner, 
& Davis (1998). 

We performed the regression analysis comparing the 
observed 7J-band light curve with the predicted il-band 
light curve in equation 4. We adopted only the data set 
where observations in V, H and K were carried out on 
the same night in order to decrease the effect of interpo- 
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lation. The parameters used for SF(7) = Sqv were So = 
0.033 mJy 2 and (3 = 1.1 for V, and S Q = 0.0077 mJy 2 
and /3=2.0 for K. Zero-correlation was assumed between 
the data in any two bands, and the standard error was 
used for the error of the fitted parameters. 

We carried out the regression analysis for a range of r 
from —10 to 60 days. The value of the multiple corre- 
lation coefficient was calculated as a function of r, and 
its peak was found to occur at r = 6.0 days with the 
X 2 -value of 12.5 for m = 26 (Figure 2). The reduced 
chi-square value becomes « 0.5, which indicates that the 
fit is good enough to represent the observations. The re- 
gression at t = 6.0 days yielded that the (H — K)-co\or 
temperature for the dust component is 1635+20K, and 
the index of power-law SED for the accretion disk com- 
ponent (/„ oc v a ) is a 
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Fig. 2. — Multiple correlation coefficient as a function of lag 
time t between variations in the H and K bands for the case of 
MCG+08-11-011. The coefficient value is peaked at r = 6.0 days. 



4. DISCUSSION 

The (H — K )-color temperature of 1635K here derived 
for MCG+08-11-011 agrees well with the sublimation 
temperature of 1500 — 1800K for graphite/silicate grains 
(Salpeter 1977; Huffman 1977). Hence, this component 
represents the innermost part of the dust torus heated to 
the sublimation temperature by the accretion disk. 

Following Winkler et al. (1992), we performed the flux 
variation gradient (FVG) analysis for the BV monitor- 
ing data (Figure 3). We converted the FVG (B - V)- 
color to the index of the power-law SED and obtained 
a=+0. 21+0.10. Combining this and the results in sec- 
tion 3 from the regression analysis, we found that the 
power-law SED with a = —0.1 ~ +0.4 prevails in the 
optical to near-infrared wavelength region (Figure 4). 

For comparison, we computed theoretical colors of the 
standard accretion disk model (Shakura & Sunyaev 1973) 
with a black hole mass of M = 10 7 M Q and an outer disk 
radius of R = !0 4 R g , where R g = GM/c 2 . Figure 4 
shows a as a function of frequency v for two cases of 
accretion rate M = 0.1M Q yr _1 (upper thick line) and 
O.OOlMgyr" 1 (lower thick line). These predicted SEDs 
have an approximate value of a = 1/3, and agree with 
our finding of a = -0.1 - +0.4 for MCG+08-11-011. 

The standard disk model has a critical difficulty for 



explaining the variability. Although many UV/optical 
monitoring observations have so far been made, all of 
them have failed to detect significant lags between UV 
and optical variations, except for NGC7469 (Peterson et 
al. 1998; Wanders et al. 1997; Collier et al. 1998) and 
Ark564 (Collier et al. 2001). The expected lag time is 
too short for the energy to flow at sound speed within the 
standard viscous disk. However, such a lag time is suffi- 
cient for the incident radiation to transfer the energy to 
any distant parts of the disk (Krolik et al. 1991). Hence, 
the irradiated accretion disk model has been invoked, 
where the central X-ray energy source illuminates a ge- 
ometrically thin disk, and UV/optical light is radiated 
at the outer part of the disk. The predicted spectrum 
of the irradiated standard disk model has an index of 
a = 1/3 at the radii where the height of the X-ray source 
is negligible compared to the radial distance from the 
center. This spectrum is exactly the same as that of the 
standard viscous disk model. Furthermore, the detected 
wavelength-dependent lag time agrees with r cx A 4 / 3 in 
the UV to optical region for NGC7469 and Ark564, which 
is the same as the prediction for the irradiated standard 
disk model (Collier et al. 1998; Collier et al. 2001). 
Consequently, we conclude that the power-law SED ob- 
tained by our analysis reinforces evidence that irradiated 
standard disks exist in the center of AGNs. 
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Fig. 3.— The B-flux to V-flux diagram for MCG+08-11-011. 
The straight line is the least-square fit to the data. 
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Fig. 4. — Values of power index a plotted against the rest-frame 
frequency of radiation. Filled circles are our results for MCG+08- 
11-011 from the data of BV , VH and VK . Vertical bar indicates 
the estimated error of a, and horizontal bar corresponds to the 
wavelength range between two bands used for calculation of a. For 
reference, open circles show the values of a for other AGNs by 
Winkler et al. (1992), and two thick lines show the predictions of 
two standard accretion disk models. See text for detail. 

Our analysis showed that the power-law SED extends 
to the near-infrared wavelength region. By the standard 
accretion disk model, the near-infrared /Cemitting re- 
gion in the disk is located at a radial distance of R « 
10 3 R g , from the center of typical Seyfert galaxies, which 
is close to the inner part of the broad line region (BLR) 
where high-ionization lines originate. Collin & Hurc 
(2001) showed that the disk would become gravitation- 
ally unstable at R > R crit w 2 x 10 4 R g (M/l0 7 M Q )-° A(i . 
Based on this estimation, the ^-emitting region in the 
disk could be gravitationally stable, and the existence 
of such a stable extended disk could distort the Kepler 
motion of BLR gas clouds. This may explain the signifi- 
cant scatter of observed virial relationships of the broad- 
emission lines, in accord with the suggestion of Peter- 
son et al. (2004) that the possible existence of some 
additional components may affect the BLR gas motion. 
Thus, near-infrared observations of AGNs bring infor- 
mation of physical coupling between the extended disk 
and the BLR. Particularly, in the case of AGNs with 
high Eddington ratio, the if-emitting region in the disk 
would be further away from the centers of AGNs, and 
the outer disk edges might be observable (Collin et al. 
2002; Kawaguchi 2003; Kawaguchi et al. 2004). 

The accretion disk component contributes to the near- 
infrared flux by 15-30% in the H band and by 15-25% in 
the K band (Figure 5). Since the observational errors are 
of order 1% only, these disk contributions are significant, 
and cannot be ignored in deriving accurate and reliable 
results from near-infrared observations of AGNs. For a 
long time, however, separation of accretion disk compo- 



nents in the near-infrared region has been considered dif- 
ficult, because the thermal emission from the dust torus 
component is the dominant source there. In this context, 
the multiple regression analysis for multicolor monitoring 
observations is indeed a promising method for extracting 
the accretion disk component from near-infrared data, 
and will be useful for investigating the physical condi- 
tions and emission mechanisms of both components of 
accretion disk and dust torus simultaneously. 
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Fig. 5. — Fractional contribution of near-infrared flux from the 
accretion disk component throughout the monitoring observations 
of MCG+08-11-011. Shown arc the values of F^ ar ' disk (t) / F x (t) 
for the H band (open squares) and for the K band (filled squares). 
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